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Abstract: We derive a sufficient condition for the existence of index- 
guided modes in a very general class of dielectric waveguides, including 
photonic-crystal fibers (arbitrary periodic claddings, such as "holey 
fibers"), anisotropic materials, and waveguides with periodicity along the 
propagation direction. This condition provides a rigorous guarantee of 
cutoff-free index-guided modes in any such structure where the core is 
formed by increasing the index of refraction (e.g. removing a hole). It 
also provides a weaker guarantee of guidance in cases where the refractive 
index is increased "on average" (precisely defined). The proof is based on a 
simple variational method, inspired by analogous proofs of localization for 
two-dimensional attractive potentials in quantum mechanics. 
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1. Introduction 

In this paper, we present rigorous sufficient conditions for the existence of index-guided modes, 
including conditions for cutoff-free modes, in a wide variety of dielectric waveguides — from 
ordinary step-index fibers [1], to photonic -crystal "holey" fibers [2-5], and even fiber-Bragg 
gratings [6] or other periodically modulated waveguides [5,7,8]. The dispersion relations of 
such waveguides must almost always be computed numerically, and so exact analytical theo- 
rems like the one derived here provide a foundation of certainty that is not available in any other 
way. A rigorous theorem allows us to give a general answer (although not a necessary condi- 
tion) for questions such as: if the waveguide core has a mixture of higher- and lower-index 
regions, how much higher-index material is enough for cutoff-free guidance; and under what 
conditions do photonic-crystal fibers, like step-index fibers, have cutoff-free guided modes? 
The theorem provides an absolute guarantee, with no calculation required, that strictly increas- 
ing the refractive index to form the waveguide (e.g. filling in a hole of a holey fiber) yields a 
cutoff-free guided mode. It also leads directly to necessary conditions for single-polarization 
fibers, the subject of another manuscript currently in prepration. Our work extends an earlier 
proof of guided modes for homogeneous-cladding, non-periodic, dielectric waveguides with 
isotropic [9] or anistropic [10] materials, and is closely related in spirit to proofs of the exis- 
tence of bound modes in two-dimensional potentials for quantum mechanics [11]. 

The most common guiding mechanism in dielectric waveguides is index guiding (or "total in- 
ternal reflection"), in which a higher-index core is surrounded by a lower-index cladding e c (e is 
the relative permittivity, the square of the refractive index in isotropic non-magnetic materials). 
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(a) Cross section of a 
waveguide (e.g. a 
conventional fiber) 
with a homogeneous 
cladding and an 
arbitrary- shape core. 



(b) Cross section of a 
photonic-crystal fiber 
with periodic cladding 
and arbitrary-shape core. 



(c) A waveguide periodic 
in the propagation (z) 
direction surrounded 
by a homogenous 
cladding. 



Fig. 1 . Schematics of various types of dielectric waveguides in which our theorem is ap- 
plicable. Light propagates in the z direction (along which the structure is either uniform 
or periodic) and is confined in the xy direction by a higher-index core compared to the 
surrounding (homogeneous or periodic cladding). 



A schematic of several such dielectric waveguides is shown in Fig.[T] In particular, we suppose 
that the waveguide is described by a dielectric function e(x,y,z) — e c (x,y,z) + Ae(x,y,z) such 
that: e, e c , and Ae are periodic in z (the propagation direction) with period a {a — » for the 
common case of a waveguide with a constant cross-section); that the cladding dielectric func- 
tion e c is periodic in xy (e.g. in a photonic-crystal fiber), with a homogeneous cladding (e.g. 
in a conventional fiber) as a special case; and the core is formed by a change Ae in some re- 
gion of the xy plane, sufficiently localized that f | 1/e — 1 /e c \ <°° (integrated over the xy plane 
and the unit cell in z). This includes a very wide variety of dielectric waveguides, from con- 
ventional fibers [Fig. [TJa)] to photonic-crystal "holey" fibers [Fig. [TJb)] to waveguides with a 
periodic "grating" along the propagation direction [Fig.[TJc)] such as fiber-Bragg gratings and 
other periodic waveguides. We exclude metallic structures (i.e, we require £ > 0) and make the 
approximation of lossless materials (real e). We allow anisotropic materials. The case of sub- 
strates (e.g. for strip waveguides in integrated optics [12-14]) is considered in Sec. [5] We also 
consider only non-magnetic materials (relative permeability jj. = 1), although a future extension 
to magnetic materials should be straightforward. Intuitively, if the refractive index is increased 
in the core, i.e. if Ae is non-negative, then we might expect to obtain exponentially localized 
index-guided modes, and this expectation is borne out by innumerable numerical calculations, 
even in complicated geometries like holey fibers [2-5]. 

However, an intuitive expectation of a guided mode is far from a rigorous guarantee, and upon 
closer inspection there arise a number of questions whose answers seem harder to guess with 
certainty. First, even if Ae is strictly non-negative, is there a guided mode at every wavelength, 
or is there the possibility of e.g. a long-wavelength cutoff (as was initially suggested in holey 
fibers [15], but was later contradicted by more careful numerical calculations [16])? Second, 
what if Ae is not strictly non-negative, i.e. the core consists of partly increased and partly 
decreased index; it is known in such cases, e.g. in "W-profile fibers" [17] that there is the 
possibility of a long-wavelength cutoff for guidance, but precisely how much decreased-index 
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Fig. 2. Example dispersion relation of a simple 2d dielectric waveguide in air (inset) for the 
TM polarization (electric field out of the plane), showing the light cone, the light line, the 
fundamental (cutoff-free) guided mode, and a higher-order guided mode with a cutoff. 

regions does one need to have such a cutoff? Third, under some circumstances it is possible to 
obtain a "single-polarization" fiber, in which the waveguide is truly single-mode (as opposed 
to two degenerate polarization modes as in a cylindrical fiber) [18-23] — our theorem can be 
extended, similar to Ref. 9, to a condition for two guided modes, and we will explore the 
consequences for single-polarization fibers in a subsequent paper. It turns out that all of these 
questions can be rigorously answered (in the sense of sufficient conditions for guidance) for the 
very general geometries considered in Fig.[T] without resorting to approximations or numerical 
computations. 

We will proceed as follows. First, in Sec. [2] we review the mechanism of index guiding, state 
our result (a sufficient condition for the existence of index-guided modes), and discuss some 
important special cases. In Sec. [3] we first prove this theorem for the simplified special case 
of a homogeneous cladding e c , where the proof is much easier to follow. Then, in Sec. [4] we 
generalize the proof to arbitrary periodic claddings, such as for holey photonic-crystal fibers 
(with some algebraic details left to the appendix). In Sec. [5] we discuss a few contexts that go 
beyond the initial assumptions of our theorem: substrates, material dispersion, and finite-size 
effects. Finally, we offer some concluding remarks in Sec.[6]discussing future directions. 

2. Statement of the theorem 

First, let us review the basic description of the eigenmodes of a dielectric waveguide [1, 5]. 
In a waveguide as defined above, the solutions of Maxwell's equations (both guided and non- 
guided) can be written in the form of eigenmodes H(x,y 1 z)e'P z ~" ))t (via Bloch's theorem thanks 
to the periodicity in z) [5], where CO is the frequency, j3 is the propagation constant, and the 
magnetic-field envelope H(x,y,z) is periodic in z with period a (or is independent of z in the 
common case of a constant cross section, a — > 0). A plot of CO versus j3 for all eigenmodes is the 
"dispersion relation" of the waveguide, one example of which is shown in Fig. [2] In the absence 



of the core (i.e. if Ae = 0), the (non-localized) modes propagating in the infinite cladding form 
the "light cone" of the structure [2-5]; and at each real j3 there is a fundamental (minimum-O)) 
space-filling mode at a frequency ft) c (/3) with a corresponding field envelope H t . [2-5]. Such a 
light cone is shown as a shaded triangular region in Fig. [2] Below the "light line" co c ([5), the 
only solutions in the cladding are evanescent modes that decay exponentially in the transverse 
directions [2-5,24]. Therefore, once the core is introduced (Ae ^ 0), any new solutions with 
CO < co c must be guided modes, since they are exponentially decaying in the cladding far from 
the core: these are the index-guided modes (if any). Such guided modes are shown as lines 
below the light cone in Fig. [2] in this case, both a lowest-lying ("fundamental") guided mode 
with no low-frequency cutoff (although it approaches the light line asymptotically as CO — ► 0) 
and a higher-order guided mode with a low-frequency cutoff are visible. Since a mode is guided 
if co < co c , we will prove the existence of a guided mode by showing that CO has an upper bound 
< CO c , using the variational (minimax) theorem for Hermitian eigenproblems [5]. 

[Modes that lie beneath the light light are not the only type of guided modes in microstruc- 
tured dielectric waveguides. While all the guided modes in a traditional, homogeneous-cladding 
fiber lie below the light line and are confined by the mechanism of index-guiding, there can 
also be bandgap-guided modes in photonic-crystal fibers [2-5]. These bandgap-guided modes 
lie above the cladding light line and are therefore not index-guided. Bandgap-guided modes 
always have a low-frequency cutoff (since in the long-wavelength limit the structure can be 
approximated by a "homogenized" effective medium that has no gap [25]). We do not con- 
sider bandgap-guided modes in this work; sufficient conditions for such modes to exist were 
considered by Ref. 26.] 

We will derive the following sufficient condition for the existence of an index-guided mode 
in a dielectric waveguide at a given j3 : a guided mode must exist whenever 

^•(e^-^D^O, (1) 

where the integral is over xy and one period in z and D c is the displacement field of the 
cladding's fundamental mode. From this, we can immediately obtain a number of useful special 
cases: 

• There must be a cutoff-free guided mode if Ae > everywhere (i.e., if we only increase 
the index to make the core). 

• For a homogeneous cladding (and isotropic media), there must be a cutoff-free guided 
mode if [( 1/e — 1 /e c ) < (similar to the earlier theorem of Ref. 9, but generalized to 
include waveguides periodic in z and/or cores Ae that do not have compact support). 

• More generally, a guided mode has no long-wavelength cutoff if eq. ([TJ is satisfied for 
the quasi-static (CO — > 0, j3 — > 0) limit of D c . 

Equation ([1} can also be extended to a sufficient condition for having two guided modes (or 
equivalently, a necessary condition for single-polarization guidance), when the cladding fun- 
damental mode is doubly degenerate. We explore this generalization, analogous to a result in 
Ref. 9 for homogeneous claddings, in another manuscript currently being prepared. 

3. Waveguides with a homogeneous cladding 

To illustrate the basic ideas of the proof in a simpler context, we will first consider the case 
of a homogeneous cladding (e 6 = constant) and isotropic, z-invariant strucures (e is a scalar 
function of x and y only). In doing so, we reproduce a result first proved (using a somewhat 
different approach) by [9] (although the latter result required Ae to have compact support, 



whereas we only require a weaker integrability condition). Our proof, which we generalize in 
the next section, is closely inspired by a proof [11] of a related result in quantum mechanics, 
the fact that any attractive potential in two dimensions localizes a bound state [11,27-30]; we 
discuss this analogy in more detail below. 

That is, we take the dielectric function e(x,y) to be of the form: 

e(x,y) = e c +Ae(x,y), (2) 

where Ae is an an arbitrary change in e that forms the core of the waveguide. For convenience, 
we define a new function A by: 

A(x,y)^e- l -e- 1 . (3) 

The only constraints we place on Ae are that £ be real and positive and that / \A\dxdy be finite, 
as discussed above. Now, we wish to show that there must always be a (cutoff-free) guided 
mode as long as Ae is "mostly positive," in the sense that: 

J A(x,y)dxdy < 0, (4) 

Since eq. Q is independent of CO or j3, the existence of guided modes will hold at all frequencies 
(cutoff -free). 

The foundation for the proof is the existence of a variational (minimax) theorem that gives 
an upper bound for the lowest eigenfrequency ©, n m. m particular, at each fi, the eigenmodes 
13.{x,y)e'P z ~ iat satisfy a Hermitian eigenproblem [5]: 

V /3 xiv /3 xH-0 i 3H=^H, (5) 

where 

+ (6) 

with eq. |5j) defining the linear operator . In addition to the eigenproblem, there is also the 
"transversality" constraint [5,24]: 

V r H = (7) 

(the absence of static magnetic charges). From the Hermitian property of 0^, the variational 
theorem immediately follows [5]: 

c 2 v r H=o JH*-Hdxdy 

That is, an upper bound for the smallest eigenvalue is obtained by plugging any "trial function" 
H(x,y), not necessarily an eigenfunction, into the right-hand-side (the "Rayleigh quotient"), as 
long as H is "transverse" [satisfies eq. ^] Q Conversely, if eq. is not satisfied, it is easy 
to make the numerator of the right-hand-side (which involved x H) zero, e.g. by setting 
H = V<p + ificpz, for any cp(x,y), so transversality of the trial function is critically important to 
obtaining a true upper bound. 

Now, we merely need to find a transverse trial function such that the variational upper bound 
is below the light line of the cladding, which will guarantee a guided fundamental mode. For a 



'Technically, we must also restrict ourselves to trial functions where the integrals in eq. jis} are defined, i.e. the trial 
functions must be in the appropriate Sobolev space //(V p x). 



homogeneous, isotropic cladding £ c , the light line is simply (0 2 /c 2 = j3 2 /e t , and so the condi- 
tion for guided modes becomes: 



e c J H* • ©pHdxdy - /r J H* • Bdxdy 

= e c J ^\\VpxH\\ 2 dxdy-l3 2 J \\H\\ 2 dxdy <0, (9) 

where in the second line we have integrated by parts. 

The problem of bound states in quantum mechanics is conceptually very similar. There, given 
a potential function V(x,y) in two dimensions with f\V\ <°°, one wishes to show that JV < 
(attractive) implies the existence of a bound state: an eigenfunction of the Schrodinger operator 
— V 2 + V with eigenvalue (energy) < 0. Again, this is a Hermitian eigenproblem and there is 
a variational theorem [31], so one merely needs to find some trial wavefunction y/ for which 
the Rayleigh quotient is negative in order to obtain a bound state. In one dimension, finding 
such a trial function is simple — for example, an exponentially decaying function e~ a l v l (or a 
Gaussian e~ ax ) will work for sufficiently small a — and the proof is sometimes assigned as an 
undergraduate homework problem [32]. In two dimensions, however, finding a trial function 
is more difficult — in fact, no function of the form f(ocr) (where r is the radius \J x 2 +y 2 ) will 
work (without more knowledge of the explicit solution for V) [11] — and the earliest proofs of 
the existence of bound modes used more complicated, non-variational methods [27,30]. How- 
ever, an appropriate trial function for a variational proof was eventually discovered [9,29], and 
later a simpler trial function e~( r+1 ) was proposed independently by Yang and de Llano [11], 

In the present electromagnetic case, we found that the following trial function, inspired by 
the quantum case above [11], works. That is, we can prove the existence of waveguided modes 
for a homogeneous cladding using the trial function, in cylindrical (r, 0) coordinates: 

H = f ycos — (ry)' sin , (10) 

where 

Y=Y(r)=e l - (r2+1)a (11) 



for some a > 0, and (ry)' is the derivative with respect to r. Clearly, H in eq. ( 10 1 reduces to an 
x-polarized plane wave propagating in the z direction as a — > (and hence y — > 1). This is a key 
property of the trial function: in the limit of no localization (a = 0, Ae = 0) it should recover a 
fundamental (lowest-o) solution of the infinite cladding. Also, by construction, it satisfies the 
transversality condition |7]i (which is why we chose this particular form). We chose y slightly 
differently from Ref. 1 1 for convenience only (to make sure it is differentiable at the origin and 
goes to 1 for a — > 0). For future reference, the first two r derivatives of y are: 

= -2ar{r* + l) a ~ l y, (12) 
f = 2a(r 2 + l) a - l Y [-1 +2ar 2 (r 2 + l)"" 1 +2(1 - a)r 2 (r 2 + l)" 1 ] , (13) 

and are plotted along with y in Fig. [3] 

What remains is, in principle, merely a matter of algebra to verify that this trial function, for 
sufficiently small a, satisfies the variational condition ([9j. In practice, some care is required in 
appropriately bounding each of the integrals and in taking the limits in the proper order, and we 
review this process below. 
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Fig. 3. Plot of y [eq. (IT}], //a [eq. (l2)], and f/a [eq. (IJ}] versus r for a = 0.1. All 
three functions go to zero for r — > o°, with no extrema other than those shown. 



We substitute the trial function (jTOj for H into the left-hand side of eq. |9}: 
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We proceed to show that the last line of the above expression is negative in the limit a — > 0, 
thus satisfying the condition for the existence of bound modes. We first examine the second 
term of eq. (14) : 



lim 



/3 2 A(r)||H! 



j3 2 / A(ry 2 r < 0. 



(15) 



The key fact here is that we are able to interchange the a — > limit and the integral in this case, 
thanks to Lebesgue's Dominated Convergence Theorem (LDCT) [33]: whenever the absolute 
value of the integrand is bounded above (for sufficiently small a) by an a-independent function 
with a finite integral, LDCT guarantees that the a — > limit can be interchanged with the 
integral. In particular, the absolute value of this integrand is bounded above by |A| multiplied 
by some constant (since |H| is bounded by a constant: |y| < 1 and |r/| is also easily seen to 
be bounded above for sufficiently small a), and |A| has a finite integral by assumption. Since 
lim a ^o |H| 2 = 1, we obtain eq. (Wb, which is negative by assumption. 

Now we must show that the remaining first term of eq. ( 14 1 goes to zero as a — > 0, completing 
our proof. This term is proportional to (g^ 1 +A), but the A terms trivially go to zero by the same 
arguments as above: A allows the limit to be interchanged with the integration by LDCT, and 
as a — + the / and y" terms go to zero. The remaining e^T 1 terms can be bounded above by a 



sequence of inequalities as follows: 
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From the first to second line, we substituted eqs. ( |T2] > and ( 13 i and simplified. From the second 
to third line, we bounded the integrand above by flipping negative terms into positive ones 
and replacing r 2 with r 2 + 1 . From the third to the fourth line, we made a change of variables 



t 2 = r z + 1 . Then, from the fourth to fifth line, we made another change of variable u = t 2a , and 
bounded the integral above by changing the lower limit from u = 1 to u = 0. The final integral 
can be performed exactly and goes to zero, completing the proof. 

4. General periodic claddings 

In the previous section we considered z-invariant waveguides with a homogeneous cladding 
and isotropic materials (for example, conventional optical fibers). We now generalize the proof 
in three ways, by allowing: 

• transverse periodicity in the cladding material (photonic-crystal fibers), 

• a core and cladding that are periodic in z with period a {a — > for the z-invariant case), 

• anisotropic e c and Ae materials (e is a 3 x 3 positive-definite Hermitian matrix). 
In particular, we consider dielectric functions of the form: 



e(x,y,z) = e c (x,y,z)+Ae(x,y,z), 



(17) 



where the cladding dielectric tensor e c (x,y,z) = £ c (x,y,z + a) is z-periodic and also periodic 
in the xy plane (with an arbitrary unit cell and lattice), and the core dielectric tensor change 
Ae(x,y,z) = Ae(x,y,z + a) is z-periodic with the same period a. Both e c and the total e must 
be positive-definite Hermitian tensors. As defined in eq. ([3]), we denote by A the change in 
the inverse dielectric tensor. Similar to the isotropic case, we require that / |A,- ; | be finite for 
integration over the xy plane and one period of z, for every tensor component Ay. We also 
require that the components of e^T 1 be bounded above. 

In the homogeneous -cladding case, any light mode that lies beneath the (linear) light line 
of the cladding is guided. We have shown that such a mode always exists, for all j3, under the 
condition of eq. Q, by showing that the variational upper bound on its frequency lies below 
the light line. In the case of a periodic cladding, the light line is the dispersion relation of the 
fundamental space-filling mode of the cladding, which corresponds to the lowest-frequency 
mode at each given propagation constant j5 [2-5]. This light "line" is, in general, no longer 
straight, and there are mechanisms for guidance that are not available in the previous case, 



such as bandgap guidance [2-5]. Bandgap-guided modes may exist above the light line and 
are, in general, not cutoff-free because the gap has a finite bandwidth. Here, we only consider 
index-guided modes, which are guided because they lie below the light line. We will follow 
the same general procedure as in the previous section to derive the sufficient condition [eq. ([TJ] 
to guarantee the existence of guided modes. The homogeneous -cladding case is then a special 
case of this more general theorem, recovering eq. Q (but generalizing it to z-periodic cores), 
where in that case the cladding fundamental mode D 6 is a constant and can be pulled out of 
the integral. The case of a z-homogeneous fiber is just the special case a — + 0, eliminating the z 
integral in eq. ([TJ. 

The proof is similar in spirit to that of the homogeneous-cladding case. At each j3, the eigen- 
modes H(x,y,z)e'^~ ,<i)r satisfy the same Hermitian eigenproblem ^ and transversality con- 
straint ([TJ as before. We have a similar variational theorem to eq. ([8} [5], except that, in the case 
of z-periodic ity, we now integrate over one period in z as well as over x and y. 

c z Vp-H=o JH* H 

As before, to prove the existence of a guided mode we will find a trial function H such that 
this upper bound, called the "Rayleigh quotient" for H, is below the light line (O t (j3) 2 /c 2 . The 
corresponding condition on H can be written [similar to eq. {9}]: 



/ 



H* ©pH- I H H 0. (19) 



We considered a variety of trial functions, inspired by the Yang and de Llano quantum 
case [11], before finding the following choice that allows us to prove the condition (19) . Sim- 
ilar to eq. ( lOl, we want a slowly decaying function proportional to y(r) = e 1_ (' 2 + 1 ) , from 



eq. (11 1, that in the a — > (weak guidance) limit approaches the cladding fundamental mode 
H c . As before, the trial function must be transverse (V^ • H = 0), which motivated us to write 
the trial function in terms of the corresponding vector potential. We denote by A c the vector 
potential corresponding to the cladding fundamental mode H e = x A c . In terms of A c and 
y, our trial function is then: 

H = Vg x (yA c ) = yH, + Vy x A,. (20) 

For convenience, we choose A e to be Bloch-periodic (like H c , since A c also satisifies a pe- 
riodic Hermitian generalized eigenproblem and hence Bloch's theorem applies)]^] In contrast, 



our previous homogeneous-cladding trial function [eq. ( 10 1] corresponds to a different gauge 
choice with an unbounded vector potential A ( = — ^ y + y/, differing from a constant vector 

potential via the gauge function \j/ = 4g sin0 +e~'P z . 



Substituting eq. (20 1 into the left hand side of our new guidance condition ( 19 1, we obtain 



five categories of terms to analyze: 

(i) terms that contain A = e — e~ l , 

(ii) terms that cancel due to the eigenequation Q, 

(iii) terms that have one first derivative of 7, 

(iv) terms that have (/) 2 , 



2 Alternatively, it is straightforward to show that the Coulomb gauge choice, Vp ■ A c = 0, gives a Bloch-periodic A c , 
by explicitly constructing the Fourier-series components of A c in terms of those of H r . 



(v) terms that have /y" or (y") 2 . 



Category (i) will give us our condition for guided modes, eq. ([TJ, while category (ii) will be 
cancelled exactly in eq. ( fL9| ). We show in the appendix that all of the terms in category (iii) 
exactly cancel one another. The terms in categories (iv) and (v) all vanish in the a — ► limit; 
we distinguish them because category (v) turns out to be easier to analyze. There are no terms 
with y" alone, as these can be integrated by parts to obtain category (iii) and (iv) terms. In the 
appendix, we provide an exhaustive listing of all the terms and how they combine as described 
above. In this section, we only outline the basic structure of this algebraic process, and explain 
why the category (iv) and (v) terms vanish as a — > 0. 
Category (i) consists only of one term: 



limy H* • (Vjb x AV^ x H) 
H 6 * • (V p x AV p x H,) 



(21) 

(V^xH c )*-A(V^xH c ) 



CO? 



- c fy D: .AD, 

From the first to the second line, we interchanged the limit with the integration, thanks to the 
LDCT condition as in Sec. [3] since the magnitudes of all of the terms in the integrand are 
bounded above by the tensor components |Ay| multiplied by some a -independent constants, 
and | Ay | has a finite integral by assumption. (In particular, the A c fundamental mode and 
its curls are bounded functions, being Bloch-periodic, and y and its first two derivatives are 
bounded for sufficiently small a.) The result is precisely the left-hand side of eq. which is 
negative by assumption. 

Next, we would like to cancel — ^ J H* • H by the eigen-equation (j^J. Thus, we examine the 
term /H* • (Vp x e^yVjS x H c ) (which comes from the term where the right-most curl falls 
on H c rather than y) below: 



H* • (yV^ x e^Wp x H c + (Vy) x e" 1 x H c ) 



co 2 r (22) 

H*-y-^H c + H*-(V Y xe- l VpxH c ) 



H* • -§- H - J H* • -f- Vy x A c + J H* • (Vy x g" 1 x H c ) 

From the second to the third lines, we used the eigenequation (|5j, and from the third to the 
fourth lines we used the definition d20| of H in terms of H c . The first term of the last line 



above cancels — -j- J H* • H in eq 

-,>2 
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The second and third terms contain two category (iii) 

terms: % / yH c . • (Vy x A c ) and ico c J yVy • E c x H*, both of which will be exactly cancelled as 
described in the appendix, as well as some category (iv) and (v) terms. 

The category (iv) integrands are all of the form (y') 2 multiplied by some bounded function 
(a product of the various Bloch-periodic fields as well as the bounded e^T 1 ). This integrand 
can then be bounded above by replacing the bounded function with the supremum B of its 
magnitude, at which point the integral is bounded above by 2nB J^(Y) 2 rdr. However, such 



integrands were among the terms we already analyzed in the homogeneous -cladding case, in 
eq. ( [To} , and we explicitly showed that such integrals go to zero as a — > 0. 

The category (v) integrands could also be explicitly shown to vanish as a — > 0, similar to 
eq. ( 16 1, but a simpler proof of the same fact can be constructed via the LDCT condition. In par- 
ticular, similar to the previous paragraph, after replacing bounded functions with their suprema 
we are left with cylindrical-coordinate integrands of the form //'r and (y") 2 r. Both of these 
integrands, however, are bounded above by an a-independent function with a finite integral, 
and hence LDCT allows us to put the a — > limit inside the integral and set the integrands to 



zero. Specifically, by inspection of eqs. ( 12 1 and ( 13 1, \H'\r < 4r 2 (l + 2 + 2)/(r 2 + 1) 2 ~ 6 and 
[ffr <4r(l+2 + 2) 2 /(r 2 + l) 2 - 6 for a < 8/4, and both of these upper bounds have finite 
integrals, if we take 8 to be some number < 1 /2, since they decay faster than 1 jr. 

In summary, we have shown that, if eq. ([TJ is satisfied, then the variational upper bound for 
our trial function [eq. ( [20| ] is below the light line, and therefore an index-guided mode is guar- 
anteed to exist. The special cases of this theorem, as discussed in the introduction, immediately 
follow. 



5. Substrates, dispersive materials, and finite-size effects 

In this section, we briefly discuss several situations that lie outside of the underlying assump- 
tions of our theorem: waveguides sitting on substrates, dispersive (©-dependent) materials, and 
finite-size claddings. 

An optical fiber is completely surrounded by a single cladding material, but the situation 
is quite different in integrated optical waveguides. There, it is common to have an asymmet- 
rical cladding, with air above the waveguide and a low-index material (e.g. oxide) below the 
waveguide, such as in strip or ridge waveguides [12-14]. In such cases, it is well known that 
the fundamental guided mode has a low-frequency cutoff even when the waveguide consists of 
strictly nonnegative Ae [12, 14]. This does not contradict our theorem because we required the 
cladding to be periodic in both transverse directions, whereas a substrate is not periodic in the 
vertical direction. 

We have also assumed non-dispersive materials in our proof. What happens when we have 
more realistic, dispersive materials? Suppose that e depends on CO but has negligible absorption 
(so that guided modes are still well-defined). For a given co, we can construct a frequency- 
independent e structure matching the actual e at that CO, and apply our theorem to determine 
whether there are guided modes at co. The simplest case is when Ae > for all CO, in which case 
we must still obtain cutoff -free guided modes. The theorem becomes more subtle to apply when 
Ae < in some regions, because not only must one perform the integral of eq. ([TJ) to determine 
the existence of guided modes, but the condition (JTJ) is for a fixed j3 while the integrand is for a 
given frequency, and the frequency of the guided mode is unknown a priori. 

Finally, any real structure has a finite cladding. Both numerically and experimentally, this 
makes it difficult to study the long-wavelength regime because the modal diameter increases 
rapidly with wavelength (i.e. the frequency approaches the light line and the transverse de- 
cay rate becomes very slow) — in fact, it seems likely that the modal diameter will increase 
exponentially with the wavelength. In quantum mechanics (scalar waves) with a potential 
well of depth V, the decay length of the bound mode increases as e c l v when V — > 0, for 
some constant C [11,27]. In electromagnetism, for the long wavelength limit, a homogenized 
effective-medium e description of the structure becomes applicable [25], and in this effective 
near-homogeneous limit the modes are described by a scalar wave equation with a "potential" 
— 2 Ae [34], and hence the quantum analysis should apply. Thus, by this informal argument, 
we would expect the modal diameter to expand proportional to <? for some constant C (where 
A = 2kc/co is the vacuum wavelength), but a more explicit proof would be desirable. 



6. Concluding remarks 

We have demonstrated sufficient conditions for the existence of cutoff-free guided modes for 
general microstmctured dielectric fibers, periodic in either or both the z direction and in the 
transverse plane. The results are a generalization of previous results on the existence of such 
modes in fibers with a homogeneous cladding index. Our theorem allows one to understand 
the guidance in many very complicated structures analytically, and enables one to rigorously 
guarantee guided modes in many structures (especially those where Ae > everywhere) by 
inspection. There remain a number of interesting questions for future study, however, some of 
which we outline below. 

Our eq. (jTJ is a sufficient condition for index-guided modes, but it is certainly not neces- 
sary in general: even when it is violated, one can have guided modes with a cutoff (as for 
W-profile fibers [17] or waveguides on substrates [12, 14]), or other types of guided modes 
(such as bandgap-guided modes [2-5]). However, these other types of guides modes in di- 
electric waveguides have a long-wavelength cutoff, so one can pose the question: is eq. ([T| a 
necessary condition for cutoff-free guided modes (where D c is given by the long-wavelength 
limit of the cladding fundamental mode) in dielectric waveguides (as opposed to TEM modes 
in metallic coaxial waveguides, which also have no cutoff [35])? Based on theoretical reasoning 
and some numerical evidence, we suspect that the answer is no, but that it may be possible to 
modify eq. ([T} to obtain a necessary condition. In particular, the variational theorem is closely 
related to first-order perturbation theory: if one has a small perturbation Ae and substitutes 
the unperturbed field into the Rayleigh quotient, the result is the first-order perturbation in the 
eigenvalue. However, when Ae is large, even if the volume of the perturbation is small, pertur- 
bation theory requires a correction due to the electric-field discontinuity at the interface [36]. 
In the long-wavelength limit, perturbation theory is corrected by computing the quasi-static po- 
larizability of the perturbation [36], and we conjecture that a similar correction to our trial field 
may allow one to derive a necessary condition for the absence of a cutoff. Equation ([TJ is still 
a sufficient condition (the variational theorem still holds even with a suboptimal trial function), 
but the preceding considerations predict that it will become farther from a necessary condition 
for the absence of a cutoff as Ae is increased, and this prediction seems to be confirmed by 
preliminary numerical experiments with W-profile fibers. 

Let us also mention five other interesting directions to pursue. First, Ref. 9 actually proved a 
somewhat stronger condition than eq. ([T]i for homogeneous claddings, in that they showed the 
existence of guided modes when the integral was < (and Ae > in some region) rather than < 
as in our condition. Although the = case seems unlikely to be experimentally or numerically 
significant, we suspect that a similar generalization should be possible for our theorem (re- 
weighting the integrand to make it negative and then taking a limit as in Ref. 9). Second, as 
discussed in Sec. [5] it would be desirable to develop a sufficient condition at a fixed co rather 
than at a fixed j3, although we are not sure whether this is possible. Third, one would like a more 
explicit confirmation of the argument, in Sec. [5] that the modal diameter should asymptotically 
increase exponentially with the square of the wavelength. Fourth, it might be interesting to 
consider the case of "Bragg fiber" geometries consisting of "periodic" sequences of concentric 
layers [37], which are not strictly periodic because the layer curvature decreases with radius. 
Finally, as we mentioned in Sec. [2] it is possible to extend the theorem to a condition for two 
guided modes in many cases where the cladding fundamental mode is doubly degenerate, and 
we are currently preparing another manuscript describing this result along with conditions for 
truly single-mode ("single-polarization") waveguides. 
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Appendix: All Rayleigh-quotient terms 

In this appendix, we provide an exhaustive listing of all the terms that appear when the trial 
function [eq. ( 20 »] is substituted into eq. ( 19 1 (the condition to be satisfied, a rearrangement of 
the Rayleigh quotient bound). Since the terms that contain A [category (i)] were already fully 
analyzed in Sec.|4](since for these terms the limits could be trivially interchanged), we consider 
only the remaining terms involving £ e (r). More specifically, the only non-trivial term to analyze 
is the A-free part of the left-most integral in eq. ( 19 1: 

H* ■ (Vyg x e c " 1 V j g x H) 



H 



H 



(23) 



^•(Vpxe-^xtVyxA,)) 



We have already seen, in eq. ( |22) i, that the first term breaks down into a term that cancels 
-§■ f H* • H in eq. ( 19 1, via the eigen-equation, and two other terms. Removing the terms can- 
celled by the eigenequation, and substituting — /^E c for e c 7 1 x H ( (Ampere's law), we have: 



H*-(VyxA c )- 



H 



ft), 



H* • x e- 1 [Vyx H, + x (Vyx A c )] 



[yH c + VyxA c ]*-(VyxA c )-^ / yVy- (E c x H*) 



.CO c 



(Vyx A c )* • (Vyx E c ) + / (yV^ x H c )* • e" 1 [Vy x H, + V^ x (Vyx A^ 
(Vy x H,)* • e" 1 [Vy x H c + V p x (Vy x A c )] 
( V p x Vy x Ac) * • e c 7 1 [ Vy x H c + V^ x ( Vy x A,)] 



ft),: 



.COc 



yH:-(VyxA c )- 



ft); 



J ||VyxA,| 



2 .®c 

— I 

c 



yVy- (E c xH;) 



c 

.COc 



.COc 



yE*-(VyxH c ) 



(VyxA c )*-(VyxE c ) + /- c 

y (Vp x E c )* ■ (Vy x A c ) + f (Vy x E c )* • (Vy x A c ) 
(VyxHcf-e-^VyxHc) 
(VyxH,)*-^ 1 [V p x (VyxA c )] +c.c. 
+ / (V J3 x(VyxA c ))*.e- 1 (V J3 x(VyxA c )) 



(24) 



Above, the first "=" step is obtained by substituting the trial function for H, integrating some 
of the Vjg x operators by parts, and distributing the derivatives of yH c by the product rule. The 
second step is obtained by using Ampere's law again, combined with integrations by parts and 
the product rule; "c.c." stands for the complex conjugate of the preceding expression. Continu- 
ing, we obtain: 

2 2 

°' K '^.(VyxA,,)-^ J ||V 7 xA,|| 2 -2/^ J yVy- 9t{E e x H c *} 
' ^ ^/ 7 H : .(V 7 xA c ) 

(25) 



c 



i-± / (V 7 xA t .)-(V 7 xE c )*+c.c. +-f / 7 H t *-(V 7 xA c ) 



J (V 7 x He)* • e- 1 (V 7 xH c )+y (V 7 x H c )* ■ e' 1 [V p x (V 7 x A c )] + c.c. 
(V /J x(V 7 xA 6 .))*- e - 1 (V /J x(V 7 xA 6 .)) 



In obtaining this expression, we have grouped terms into complex-conjugate pairs and used 

2 

Faraday's law to replace x E c with z'^rH c . At this point, we have two % / 7 H*. • ( V 7 x A,.) 
terms that exactly cancel. All of the remaining terms, except for —i^r J 7 V 7 - 3i {E c x H*}, are 
multiples of two first or higher derivatives of 7 , corresponding to category (iv) and (v) terms, 
which we proved to vanish in Sec.|4] 

The only remaining term is the E c x H* term, in category (iii). This term is identically zero 
(for any a > 0) because it is purely imaginary, whereas all of the other terms are purely real 
and the overall expression must be real. More explicitly: 

-2i^ / 7 V 7 .5R{E 6 .xH:.} 



c 

= -2/^|v 7 2 -9l{E c xH*} (26) 
= -/y/v-( 7 2 ^{E t .xH:.})+/^| 7 2 V.(5R{E c .xH:}) 

The first term of the last line is zero by the divergence theorem (transforming it into a surface 
integral at infinity), since y — > at infinity. For the second term, the integrand is the divergence 
of the time-average Poynting vector 3i {E c x H*}, which equals the time-average rate of change 
of the energy density [34], which is identically zero for any lossless eigenmode (such as the 
cladding fundamental mode). 



